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Abstract. We study the LHC phenomenology of two mass-degenerate heavy gauge 
bosons with the same quantum numbers as the Z and 7. We give a leading-order 
estimate of the number of events expected in Drell-Yan processes in terms of the 
i parameters of the model (mass and coupling) and of the LHC machine specifications 

Oh| (integrated luminosity and energy). We consider the feasibility of measuring a forward- 

^ 1' backward asymmetry for various choices of the parameters and estimate the potential 

^ . reach. We comment on how the results may affect future collider design and apply our 

results to a specific model of electro-weak symmetry breaking by way of example. 



> 

CN . 1. Introduction 

in 

I The main purpose of the LHC is to discover the origin of electro-weak (EW) symmetry 
O ' breaking and any associated new physics at the TeV scale. 



One clean signal for the LHC would be the discovery of new, higher-mass copies of 
the standard-model (SM) gauge fields. For such a scenario, the LHC reach in observable 
masses and couplings is of interest. There are theoretical motivations for such searches. 
^ I Many weakly-coupled extensions of the SM gauge group yield new heavy vector fields [1]. 

Such new fields also appear in little Higgs models [2]. Even strongly-coupled extensions, 
such as technicolor (TC), predict the existence of towers of heavy spin-1 fields [3]. 

In TC models EW symmetry is broken by a condensate, along the lines of chiral 
symmetry breaking in QCD. Thus these models do not contain a light elementary Higgs 
particle |. As a result the most promising test of these models at the LHC is a search 
for higher-mass spin-1 resonances. In particular, a copy of the full SU (2) x U{1) realized 
as heavy gauge bosons would be distinctive. If the new bosons are nearly degenerate in 
mass then the model will satisfy the precision EW constraints [5, 6], without suppressing 
the coupling to SM fermions. This situation would imply that the new boson resonances 
overlap at the LHC. 

Once the LHC reaches its design specification, proton-proton collisions will scan 
over the energy range current data allow higher resonances to occur at. In generic 

I In some strongly coupled models, with approximate conformal symmetry, a light dilaton might be 
present with interactions very similar to the Higgs of the SM [4] . We do not consider this possibility. 
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models this is in the 1 TeV region, but in many concrete reahzations the bound is above 
the 2 - 3 TeV range [7, 8]. 

In this paper, we consider the scenario in which a complete set of four spin-1 copies 
of the SM gauge bosons with heavy degenerate masses are produced in pp collisions. 
We concentrate on detection of the neutral heavy-Z and heavy-7 through their decay 
into muons. The phenomenology of the degenerate scenario has been largely ignored to 
date in the literature. This is in contrast to the single Z' situation, see for example [9] 
and references therein. For simplicity, we model the decay width of these particles by 
assuming it to be dominated by the decay into SM fermions. This avoids complications 
(and model-dependencies) arising when including the decay into light scalars and gauge 
bosons. If the underlying dynamics is based on an SU{Nc) gauge theory with large 
values of Nc then the boson to boson decays are suppressed by the Nc counting. 

At LHC energies the SM background is mostly hadronic, so by looking at the neutral 
sector of higher resonances we exploit the very low background for muon and electron 
pair production. Within these assumptions we can study the event rate to leading order. 
Our aim is to identify what portion of parameter space will be testable at the LHC as 
a function of integrated luminosity and energy. 

In the SM, a forward-backward asymmetry arises from interference terms between Z 
and 7 exchange diagrams. It is a sizable effect, even in the TeV energy range. Measuring 
such an effect at a hadron collider is onerous but a method and its use at the LHC was 
surveyed for instance by Dittmar [10]. It was shown that a deviation of the measured 
asymmetry from the SM prediction is a possible signature of the presence of a new, very 
broad resonance with the quantum numbers of a Z'. 

In contrast we will study a relativity narrow resonance. In our case, near the pole 
of some new resonance well into the TeV region the number of events is well in excess 
of the SM. So much so that the study of the di-lepton invariant-mass distribution will 
contain a clear peak. This will be a suitable discovery signature. Furthermore for all 
practical purposes the SM contribution to this peak can be ignored. In this case, one 
expects no significant asymmetry to be measurable at the peak from the SM. 

If such a narrow peak is discovered in the data, but it results from the existence 
of a 7', in addition to the Z', then the interference between these two new states yields 
a large asymmetry. Measuring this asymmetry is hence a clean way to identify the 
presence of two distinct particles in the peak region, with masses so close together that 
they cannot be resolved in the invariant- mass distributions [11]. 

The organization of our paper is the following. We begin by introducing a model 
that encompasses the phenomenology of many specific models of EWSB. We then 
summarize some known results about decay widths and cross-sections. Then we proceed 
to study the Drell-Yan production of spin-1 resonances at the LHC. To illustrate the 
dependence of the event rate on model parameters we make a series of approximations to 
give a simple formula for the total number of events near a peak. We provide estimates 
of the discovery reach for two collider energies, 10 and 24 TeV. 

We then derive a simple formula for computing the reach of forward-backward 
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asymmetry measurement at a collider. At a 10 TeV machine the ability to search for 
a forward-backward asymmetry, of the type wc arc interested in, is limited to already 
excluded parameter space. Therefore we present data at 14 TeV (which is close to the 
border of feasibility) and 24 TeV. To illustrate the results we finish by concentrating 
on a particular model. Individual models give relations between the free parameters 
in our model-independent analysis and restrict the parameter space. This allows us to 
evaluate a typical reach scenario. 

A few words of explanation about these choices are needed. The 10 TeV centre of 
mass energy analysis is representative of the LHC experimental program. As we will see, 
large integrated luminosities are needed for discovery and hence comparing to the 7 TeV 
- with low luminosity - initial phase of LHC running is not very useful for the cases we 
study. By 10 TeV energy, however, the reach is quite substantial. This is the reason for 
our conservative choice of energy. The 24 TeV energy analysis is somewhat speculative, 
but provides a way to gauge how an energy upgrade compares to a luminosity upgrade 
of the apparatus, in reference to the specific process we study. 

2. Set-up and Rationale 

This section reviews cross-section and asymmetry calculations as well as defining the 
notations used throughout the paper. Most of the material in this section is well-known, 
but it is useful to remind the reader about it, and to fix the notation used in later 
sections. We model the new neutral resonances by adding an extra set of SU{2) x ^7(1) 
gauge fields into the SM electro-weak Lagrangian. We assume that some mechanism 
spontaneously breaks the SU{2Y x C/(l)^ gauge symmetry to the SM SU{2)l x U{1)y 
above the TeV scale, and that in the process a copy of the W, Z and 7 acquire a common 
mass m. In general, the latter is not rigorously correct: the new heavy bosons would 
acquire different masses in any realization. However, the idea that the mass splitting 
is parametrically small is suggested by the precision electro-weak parameters. We will 
provide a simple illustration of this in the next sub-section. To parametrize the couplings 
of the new heavy gauge bosons to the SM currents we introduce a parameter R. This is 
defined so that if the coupling of a given SM boson to a given SM current is qsm then 
the corresponding new heavy boson has a coupling of y/RgsM to the same current. For 
example if the SM Z coupling to an electron current is gz then the corresponding new 
boson, the Z', would have a coupling of y/Rgz to the electron current. The common 
mass m and the universal ratio it! of the couplings are the only new parameters in the 
model. 

It must be noted that in the paper we focus our attention on the neutral vectors 
only, setting aside the charged W. This choice is due to the fact that for the processes 
we consider, the kinematics of the leptonic final state can be reconstructed precisely 
when the intermediate state is neutral, while this is difficult for charged states, in which 
case one has either to deal with missing energy, or to consider final states with hadronic 
jets. 
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2.1. Mass Degeneracy 

As previously noted, in this work we assume that the new heavy gauge fields are 
degenerate in mass. As an example, illustrating the rationale behind this assumption, 
we examine the precision parameter S [5, 6]. Recall that on general grounds one can 
write the precision parameter S as 

with g the SU{2)l coupling in the SM, / and m being the decay constant and mass 
of the vector (v) and axial (a) mesons. In models where vector dominance holds the 
lightest among the states dominate the sum. Denoting them states as the techni-p 
meson and the techni-oi meson we truncate the sum to these two states. 

Under this truncation, we write the parameters of these two states in terms of a 
constant ai and a mass splitting Am^ , 

fl -fl + aA,^m\ (2) 
^ rUp 

ml^ ^ml + Am^ (3) 
Working with these definitions to leading order in Am'^/m'i we see that. 



p 



1 P Am^ 

4 mj^ 

Experimentally, S is found to be 0(10"^) [6]. This requires to implement some 
mechanism to obtain an acceptable value from eq.(4). Choosing ai — 1 would require 
some sort of fine-tuning, hence we assume that ai be some generic 0{1) number. Another 
possibility, that fp/rn^p^ 1, is realized in QCD. However, this ratio scales as Nc in the 
large- A^c limit (for A^^ the number of colors of the confining gauge theory). This makes 
it an unsuitable mechanism in a large- A^c theory. Also, allowing for fp/m'^p ^ 1 to be 
present would mean that the (techni-)rho mesons coupling to the fermions is heavily 
suppressed (the parameter R is very small), and its main decay is into longitudinally 
polarized standard-model vector mesons. This is in analogy to what happens in QCD 
where the rho meson decays into pions. Finally, in this case the rho mesons are actually 
strongly coupled, in the sense that their width is large. For all of these reasons we prefer 
not to use this suppression mechanism. Another possibility to satisfy the electro- weak 
precision constraints is to assume that Am^/m^ <^ 1, small enough that S < 10~^. 
Interestingly, this observation extends beyond the vector-dominance approximation, at 
least at large- Ac, in the sense that S is suppressed if the mass splitting between new 
axial and vector meson states is small. This is the regime which inspires this present 
work. 



2.2. Parton-level processes 



In this section we begin the computation of the relevant cross-section, following standard 
textbook procedures. Initially this is done for a lepton machine then the results are 
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extended to a hadron machine. This choice is made to highhght the physics of the 
model considered. 

We find it useful to introduce the function Py, 

Py — {s — my + imv^v)~^ : (5) 
where V§ is the total energy, in the centre of mass frame, for the collision whose cross- 
section we compute, my is the mass of the particle V and Fy is the decay width of 
the particle V ior V = 7, Z, 7', Z'. This function enters the propagator of the boson V. 
Note that this requires computing the decay widths for the two new resonances. The 
new heavy gauge bosons can also decay into SM gauge boson pairs (e.g Z' — )■ H^"*"!^"). 
However the rate of such processes is very model-dependent, and at the same time 
assuming it to dominate the decay width might imply that the new gauge bosons be 
strongly coupled, as anticipated. In this study we assume such channels are negligible. 
This is equivalent to setting the 3-point and 4-point meson vertices to zero. This 
assumption is justified by the l/N^ counting if the underlying dynamics is similar to 
that of a large- A^c gauge theory and the heavy spin-1 states are mesons of the strongly 
coupled sector, in analogy with the p and ai of a (large- iVc) QCD-like theory. 

Within these assumptions, the 7' and Z' have decay widths proportional to those of 
the 7 and Z in the SM. The constant of proportionality is the ratio between the SM and 
new gauge boson couplings. This is the parameter R introduced earlier on. Denoting 
the SU{2) coupling of the standard model by g, the Weinberg angle by 9yj, the fine 
structure constant by a (and the new heavy boson mass by m) the decay widths into 
SM fermions of the new states are, 
8 

Fy = -Ram, (6) 
^^'"^ ^-2sin^e^ + ^sin^e^V (7) 



4% cos^ Oyj \ 3 
where we summed over all the SM fermions, including the top, neglecting the fermion 
masses. Typical values of the new heavy masses are 0(1) TeV, thus the error introduced 
by neglecting fermion masses in the phase-space factors is negligible. 

As an illustration of the model, we begin by removing the complications of pp 
collisions and take the case of e"'"e~ A*^A*~- There are four possible diagrams. In 
addition to the SM channels of 7 and Z exchange there is also 7' and Z' exchange. 
Following standard notation, introduce cy and ca for the vector and axial couplings 
defined as 

cy = -2gsin2^^, (8) 

CA = T\ (9) 

where Q denotes the electric charge and the third component of weak isospin for the 
relevant fermion. We denote 

n^(s) =47rQ;gig/(P^ + i?Py), (10) 
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for Qi and Qf the charges of the initial and final fermions in the Feynman diagram. 
Using this notation the angular 2 — )• 2 partonic cross-section can be written out as, 



|a(1 + cos^^) + Scos^}, (12) 



dn 647r2 

for 9 the angle of inclination with respect to the beam and Cl the solid angle of the 
scattering process. The A and B are defined as, 

A = |n^|' + 24U{UyU*z) + (4 + clf\Ilz\^ (13) 
B = 4clUiUjU*z) + 844|nz|' (14) 

where 3? denotes the real part. 

Similarly we can write the forward-backward asymmetry Afb ^ 

( Iq " /°i ) ^'^^ da / d cos 9 

Afb = )- -. (15) 

( /o + /-I ) dcos9 da/dcos9 

We obtain, 

= If. (16) 

It is useful to combine the above expressions and rewrite the cross-section in terms 
of an angular function g{9): 

''\iO), (17) 



dn 1927r2 
with 

g{9)^ -{1 + cos'^ 9) + Afb cos 9. (18) 
8 

This form of the partonic cross-section, in terms of g{9), is the most useful form for our 
discussion of asymmetry. 

Above an energy of 500 GeV the SM contribution becomes neghgible for our 
purposes, as anticipated in the introduction. Events near the resonance of the new 
heavy states, s ~ m^, are dominated by the new states as long as mFy/ <^ s, since in 
this case 

\Pv\ - r\ (19) 
\Pv'\ ~ (mry/)-\ (20) 

Therefore unless R is very small the new resonances will dominate. The definitions of 
Ily given can now be approximated by dropping the Py terms. 

Hy OC Py. (21) 

We use this approximation to re-write the asymmetry near the new resonance (s ~ m^) 



as 



3 clTyTz> + 24cm, 

2 r|, + 2clTyTz' + (4 + 4)'rv ■ 
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This approximation is used in computing the numerical values of Aps that we work 
with later. This method avoids the complications introduced by accounting for the 
weak dependence oi ApB on s . 

Now we quote the correct formulae in the case of a hadron collider. The formulae 
for asymmetry and cross-section at partonic level can be amended by adding a factor for 
the number of colors N^, and changing to quark couplings. Denoting the quark/lepton 
couplings with a. q/l superscript, A is modified to 

A,^i = |n^|' + 244sR(n^n|) (23) 

+ {4' + c\'){c^y' + C^/)\Uz\'. 

Including the appropriate color factor, the partonic cross-section for qq l~^l~ is 
1 s 



where the angular dependence has been integrated over. Adjusting B to account for 
quarks gives 

Bg^i = Acy^^iU^U^z) + Sclc^yc'^^lUzf. (25) 
The asymmetry at the partonic level is hence 

In summary, our study of event rates focuses around the partonic cross-section 
eq.(24). Our asymmetry calculations makes use of the angular partonic cross-section 

d'^a SsAg^i ^ 

The forward-backward asymmetry required in g{6) is computed from eq.(26) under the 
approximation of no SM background. Explicitly this means that 

2r|, + 244ryr^, + (4V4^)(ci.Vc:,')r2; 

2.3. Drell-Yan processes 

We will follow the standard recipe to model processes at a hadron collider. We consider 
the incident protons to be composed of partons holding fraction ^ of the total proton 
momentum p. Numerous sets of parton distribution functions (PDF) are available 
through Fortran code [14]. For this paper we choose the Fermi02 PDF set [15]. For 
practical purposes, one can treat the proton and its constituent partons as massless. 
The two beams of protons each have energy \fs/2. The energy of a parton-parton 
collision s is given, in terms of the two fractions and ^2, by 

s = U2S. (29) 
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If, for a parton x, the PDF is denoted by fx{C) we can follow convention and write the 
pp ^ l^l~ cross-section as, 



N 



g 

To fix the total energy ^/s of the pp collision, we impose the condition of Eq. (29). This 
means we may define the parton luminosity JSf^ to be 

Jo 

and write the full cross-section summed over all quarks as, 

^-E^^r (32) 

Finally we can write the total number of events N in some range of s at a machine with 
integrated luminosity L, 

• I ^ds. (33) 
J as 

We allow ^/s to vary between 10 and 28 TeV. The coupling ratio R will be allowed 
to vary between and 1, where the gauge boson coupling is certainly perturbative. 
Then the mass range will be taken as lying between 1 and 4 TeV. Analysis will also 
require the assumption that decay widths are small compared to the resonance mass. 
Inspection of the decay widths, under the limits on R and m, shows this condition is 
satisfied. 

Asymmetry measurements at a hadron collider are made more difficult (in 
comparison with a leptonic machine) by the unknown directions of the initial quark and 
anti-quark. However, at high rapidity there is a much greater chance that the quark 
carries the bulk of the momentum. Thus if data near a resonance in the invariant-mass 
distribution are cut so that only events with a rapidity Y greater than some value Yc 
{\Y\ > Yc) are used, then there is a corresponding degree of certainty regarding the 
quark directions. This method is used and explained more fully by Dittmar [10]. In 
summary, the magnitude of the rapidity controls the probability of misidentification 
of the quark/anti-quark directions, while the sign of the rapidity represents whether 
the quark was from the left or right in the collider. This methodology of inserting cuts, 
Lorentz boosting with the rapidity modulus to the muon centre of mass frame and finally 
using the sign of Y to define the angle consistently allows one to infer an asymmetry. 
Such an asymmetry is the remnant of the partonic asymmetry. 
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3. Model Independent Study 

3.1. Event Rate 

In this section we introduce a set of approximations that are useful in order to study the 
dependence of observable quantities (such as number of expected events and forward- 
backward asymmetry) on the parameters of the model and on the machine specifications. 
The numerical part of the analysis will, however, use the complete equations and PDF 
distributions without approximations. Prom the approximations made we gain insight 
into how event rates behave as a function of model parameters. 

Using standard electro-weak relations the coupling g can be written in terms of a 
and sin^ 9^^ ~ 0.22. Hence, near the new resonance: 

^~-^g, (^l-2sin2^^ + ^sin^^^^ ~ -Q, . (34) 

Using Eqs. (23)-(24) we can comment on the relative strength of the 7' and Z' 
contributions to the total cross-section. Because -C 1, the cross term in Eq. (24) 
is suppressed. Also, in the last term of this equation Cy^ <^ cJ^ ~ 1/4, and hence the 
contribution of the 7' is a factor of ~ IQQ^ larger than the contribution coming from 
the Z'. All of this, together with the fact that the PDF for the up-quark and down- 
quarks are the most important, and that the coupling of the up quarks to 7' is four times 
stronger than the coupling of down quarks, means that the most significant contribution 
to the total event rate is given by the subprocess liiZ — >■ 7' — >■ jjf^iJr . The next largest 
contribution to the event rate is from the cross term between the new photon and Z- 
like states (which is approximately half the value of the leading contribution) and is 
responsible for the asymmetry we discuss later. 

Focusing on the ti-quark PDF, we may attempt to approximate it with some form 
of expansion. The only relevant dimcnsionless combinations are s/s and /s. As we 
are looking at the peak we can replace the square mass with s so that only one relevant 
dimcnsionless combination remains. From the definition of =Sf the overall dimension is 

If the mass of the resonance were to increase one would expect the event rate to 
decrease. It is also reasonable to think increasing the machine energy s will increase the 
event rate. We use an approximate form of. 



and consider various n so as to obey the above criteria and then numerically fit the 
constant of proportionality. The minimal n would be —2 to fulfill our expectations for 
the behavior of =Sf. Viewing our expression as the leading term in some expansion 
a polynomial can be added to increase accuracy. We find that an acceptable 
approximation is. 




(35) 




1 



(36) 



Mass-degenerate Heavy Vector Mesons at Hadron Colliders 



10 



provided y/1 G [1,4] TeV. We then go through a numerical fit of the constants Oi and 
a2- In the range ^/s G [10, 28] TcV just those two terms are needed. 
We are interested in the event rate N in the peak region 

N (X [y^ ^ds, 
J ^ds 

~ j ^u^uds. (37) 

Approximating with its value at the pole (which is allowed since the resonances are 
very narrow), and using 



S Tfl 

-ds ~ 27r— . (38) 



{s — m'^y + m^Fy Fy 



yields 



/ 



^u^ds ^^(m')R'2n^ (39) 
as i ji 

oc R^uirri^) , (40) 

where the SM background from 7 and Z has been neglected. 

Finally, the contribution of up-quarks to the total rate is, 

sLR 1 

" ■ 1 + aimVs^ + asmVs^ ' ^ ^ 

up to a proportionality constant that contains the electro- weak couplings and the 
normalization of the PDF distribution. 

To correct for the neglected contributions of dd events and Z' exchange we can 
consider the ratio of the found rate Ny^ to a more complete expression involving uu and 
dd events with 7' and Z' exchange. For notational simplicity let us absorb the coupling 
coefficients into the definition of Py briefly and write the ratio 77 as, 

T] s\Py + Pz'?ds + ^df s\Py + Pz'?ds ^ ' 

An expression similar to Eq. (36) describes, with different choices of the coefficients, 
^d- Then we find that a very good approximation is given by 

1 l + 6i(s/s)2 + 62(s/s)S, 

for some h-i and q that arc combinations of the up and down coefficients. 
The total event rate can be written, up to an overall constant, as 

N ~ N^T]. (44) 

We fix the parameters, ai, 02, ^2; ci, C2 plus an overall normalization with the direct 
numerical computation using the PDFs. We estimate the error hence introduced 
between the estimated event rate and the direct computation to be about 5%. 

Discovery of the new resonance requires detecting a minimum number of events N 
at a machine integrated luminosity L. The approximate event rate can be solved for 
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Figure 1. The ratio R between the couphngs of the new and SM boson to SM currents, 
such that the number of events N divided by the integrated luminosity L is fixed, and 
plotted as a function of the heavy boson mass at a collision energy of y/s ~ 10 TeV. 
The leftmost line is for N/L = 5 fb then reading from left to right the lines are, 
N/L = 3, 1, 1/2 fb, finishing on the rightmost line of N/L = 1/4 fb. The events N are 
for the production of muon pairs only. 

R, given m at various machine energies and integrated luminosities, hence producing 
Figs. 1 and 2. These two plots show how the experimental reach changes when improving 
the energy of a hadron collider, in comparison with increasing only its luminosity. We 
conclude this subsection with a couple of examples, to help the reader in interpreting 
the plots. 

Suppose one considers = 10 events to be necessary for discovery of a new 
resonance. Setting = 10, if one has data corresponding to a 40 fb~^ integrated 
luminosity, then the line N/L = 1/4 fb indicates the minimum coupling R required to 
see such an effect. In the case of a 10 TeV machine, for a new resonance with mass 
m = 3 TeV, then the minimum R is around 0.4. Therefore any new boson with coupling 
ratio R > 0.4 and m < 3 TeV will produce at least 10 events in the data set obtained 
with 40 fb~^ integrated luminosity. 

Changing the machine energy from 10 TeV to 24 TeV means moving from Fig. 1 to 
Fig. 2. At 10 TeV, by improving the integrated luminosity by a factor of 20, one may 
move from a reach in masses below 2.2 TeV to roughly 3.3 TeV, provided R = 1 (see 
the leftmost and rightmost lines in Fig. 1). Whilst improving the collision energy takes 
the reach to well above 4.0 TeV, even with the no improve in integrated luminosity, for 
the same R = 1 (see the leftmost curve in Fig. 2). 

3.2. Forward-Backward Asymmetry 

In the present study we limit ourselves to theoretical considerations only, so that we 
may ascertain how feasible asymmetry measurements are from a theoretical viewpoint. 
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Figure 2. The ratio R between the couphngs of the new and SM boson to SM currents, 
such that the number of events N divided by the integrated luminosity L is fixed, and 
plotted as a function of the heavy boson mass at a collision energy of y/s ~ 24 TeV. 
The uppermost line is for N/L = 5 fb then reading from top to bottom the lines are, 
N/L = 3, 1, 1/2 fb, finishing on the lowest line of N/L = 1/4 fb. The events N are for 
the production of muon pairs only. 



In fact one may argue that, due to large PDF uncertainties, a full scale simulation of 
the experiments is overkill before LHC running. It would be interesting to see a detailed 
simulation once LHC data becomes available, and new data about the PDFs is available. 

We restrict our attention to two main sources of error in the asymmetry extraction 
algorithm: the potential for misidentification of the quark direction and finite number 
of events passing the rapidity cut Y > Yc. We set aside experimental errors, coming for 
instance from the measurement of the energy of the final leptons. 

Consider the angular distribution in Eq. (18) as a probability distribution function. 
This distribution has a mean defined as 

IT 



{6) = J eg{e) sinOde = -(2 - A 
and a variance defined by 

al= [ (e^ - {df)g{d)sinede. 



FB 



(45) 



(46) 



Interestingly, ApB turns out to be just proportional to {9). It is possible to estimate the 
deviation of the mean (9) from its true value as cTp/ yjn — 1 where n is the number of 
events used to study the angular distribution. Hence we estimate the difference between 
the true and measured asymmetry due to statistics to be, 

Aa 

AA^^* = " . (47) 

7Ty/n — 1 

Here n should be thought of as representing the number of events remaining after a cut 
in rapidity is applied to the total number of events A^. 
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Secondly there is a probability of quark misidentification introduced by the 
methodology. The assumption that the initial quark direction be always the one of the 
out-going 7' or Z', implies a possible misidentification 9 n — 9. Let the probability 
of misidentification be p(Yc), then in a sample of n events np{Yc) are misidentified. 
Let Up {^3) events travel in the forward (backward) direction in the collider. Then 
n — Up + riB and the forward-backward asymmetry is defined as {rip — tib) /n in this 
notation. The worst case scenario would be that all np{Yc) events were placed va. rip 
but should have been placed in or vice versa. In this case the difference in true and 
measured asymmetry would be, 

AA^/^= ±2p(Fe) (48) 

Adding the misidentification and the statistical together we can define a 'total 
misidentification' of, 

AA^B = ^^ + 2p(i;), (49) 
TTyn — 1 

We conservatively add the two hnearly. 

Now we define P{Y<^ to be the fraction of events with \Y\ > Yc, i.e. the probability 
an event survives the chosen rapidity cut. Then we can relate the number of events in 
the peak, N , to the number used for the asymmetry extraction, n, through P{Yc) by, 

n = P{Y^)N. (50) 

Suppose we are interested in measuring an asymmetry at some level of significance. To 
do this introduce S and require that /S.ApB < S . We must impose S > 2p{Yc) so that 
the chosen cut value is capable of producing data of sufficiently low misidentification. 
This implies that we require, 

n > (——^^IJ——] + 1 (51) 



7r(<r - 2p{Y,) 



or equivalently. 



N > 



7r(<f - 2p{Y,) 



+ 1 



(52) 



total events. 

For the model considered, at the partonic level the asymmetry is ApB ^ —0.5 in 
the case of uu — > while for dd events the asymmetry would be ApB — —0.15. 

(For comparison, SM exchange plus only one new resonance the asymmetry vanishes up 
to 3 decimal places.) The measured asymmetry will then be some average of the two. 
Using what we showed before, namely that the uu events give a larger contribution to 
the total rate, we estimate the combined asymmetry to be ApB ~ —0.42. Therefore 
seeing an asymmetry at a collider would indicate that there are two particles in the new 
neutral resonance. 

To measure such an effect at a collider one would require a 'total misidentification' 
< O.IA in order to conclude, at the 3a level, an asymmetry is present. Consider the 
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case that all events happen to fall precisely at the rapidity cut Y^. Then the probability 
p{Yc) is calculable, yielding a conservative estimate: with real data having \Y\ > Y^ one 
can in principle use a larger value of p(Vc)- The fraction P{Yc) can also be computed. 

Optimal choice of Y^ for a minimal can be achieved by writing a simple 
minimization algorithm. A high Yc gives high quality data (i.e. contaminated by 
small misidentification systematics) with low piY^) but few events will pass such the 
cut, reducing also P{Yc). The converse will be true for a low cut Yc- We used an 
optimization algorithm finding a Yc such that the required number of events has been 
minimized. 

We compute the minimum number of events required for a 3cr uncertainty at given 
mass. Using the approximate event rate given by eqs. (43) and (44) we convert the 
number of events into a statement about R. The results of this process are shown in 
Figs. 3 and 4. Figure 3 shows that large increases in integrated luminosity lead to 

1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Mass (TeV) 

Figure 3. The minimum ratio R between couplings of the new and SM boson to 
SM currents, such that at the Scr level (see text for definition) a forward-backward 
asymmetry of -0.42 is detectable, plotted for various integrated luminosities at a 
collision energy of ^/s =14 TeV and as a function of the heavy boson mass. The 
leftmost line is for an integrated luminosity of 1 fb^^, the middle line for 10 fb^^ and 
the rightmost 100 fb"\ 

modest gains in reach, compared to the gain obtained by going to larger energies, as 
shown in Fig. 4. One can also see from fig. 3 that in order to test the interesting range 
m > 2 TeV, one requires a very high integrated luminosity of 100 fb~^. At a 10 TeV 
machine it is not possible to measure an asymmetry of a neutral heavy gauge boson in 
this range of most interest, and hence we omit showing the related plot. 
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4. PDF Uncertainties 

So far the uncertainty on the PDFs themselves has been completely ignored. The 
Fortran module supplied by [14] has the capability to compute PDF uncertainties. As 
the computational cost to produce uncertainties with the package is very large, we use an 
approximate expression which captures the essential features of the PDF uncertainties. 
Namely, the quark PDFs are well known up to fractions of around 0.5, then they become 
poorly understood while the anti-quark PDFs are poorly understood across all of the 
fraction range. Two piecewise functions, one for the u and one for the u, approximate 
this behavior 

An(0 _ f 0.05 
u{0 ~ I 0.8^ -0.35 

AujQ _ f 2.25^ 
u{0 ~ [25^-18.2 

These uncertainties are approximations to the results given by using the package [14] 
for our chosen PDF set [15] and the results broadly agree with those published by the 
CTEQ collaboration for their PDF sets [16]. Even if an event is not identified as being 
uu or dd in our analysis, the uu dominates events and hence the uncertainty indicated 
will offer a good guide as to the situation. 

As the asymmetry extraction methodology involves cutting the data for large Y, it 
probes the region with a smaller anti-quark fraction and a larger quark fraction. If the 
quark fraction gets above 0.5 the uncertainty grows rapidly. At a 10 TeV collider this 
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Figure 4. The minimum ratio R between couplings of the new and SM boson to 
SM currents, such that at the 3cr level (see text for definition) a forward-backward 
asymmetry of -0.42 is detectable, plotted for various integrated luminosities at a 
collision energy of -^i =24 TeV and as a function of the heavy boson mass. The 
leftmost line is for an integrated luminosity of 1 fb~^, the middle line for 10 fb^^ and 
the rightmost 100 fb"\ 



ife<0.5 

(53) 

otherwise 
if e < 0.8 

(54) 

otherwise 
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effect will come in faster as high fractions are required to make large invariant masses 
while the effect will be slower to appear in a 24 TeV collider. 

In general, most events in a peak will have low rapidity. As a result we expect far 
lower uncertainties for event rate measurements compared to asymmetry measurements. 

5. An Example Model 

In order to illustrate the results of the model-independent study it is useful to apply 
them to a specific model. For example in the AdS/TC model in [17] a 5D model is 
formulated in AdS space and bounds on the {R, m) parameter space are discussed. 

In this work, several qualitatively different regions of the (i?, m) parameter space 
are identified. Firstly there is a region where the dominant decay mode of the techni- 
p particles is into SM gauge bosons, and not SM fermions. This region is marked as 
invisible in figs. 5 & 6 because the branching ratio into fermions would strongly suppress 
the Drell-Yan process of interest here. There is then a second region already excluded 
by precision measurements. Finally the constraint R < 0.65 comes from specific details 
of the model. 

Putting this information into a plot and overlaying the search capabilities of a 
collider produces Figs. 5 and 6. At the very high integrated luminosity of Fig. 5 most 
of the parameter space will be tested for the existence of new resonances. However, 
there is no reach in terms of asymmetry. Thus it will be hard to conclude whether 
there is more than one particle responsible for the peak. In contrast, with the higher 
collision energy of Fig. 6, a very large fraction of the parameter space presents a testable 
asymmetry. 

Both Figs. 5 & 6 include a la PDF uncertainty marked by thin red lines and a 
shading between. The lower of these two fines falls below the axis and is barely visible. 
It is clear that the event rate uncertainty is much smaller in the 24 TcV case with the 
cone barely visible on the plot. As for the asymmetry, the uncertainty at 14 TcV grows 
much faster for those masses where the 14 and 24 TeV lines are both able to reach. 
In the case of the 14 TeV curve, the uncertainty in coupling is 10 times the coupling 
estimate by 2 TeV. In comparison the coupling uncertainty will not reach a factor of 10 
until 3.5 TeV in the 24 TeV machine. 

For the representative point m = 2.6 TeV and it! = 0.4 we simulate the asymmetry 
and see the rapidity cut in action. For simplicity, we assume that only uu events 
contribute at the partonic level. As explained, this is the dominant subprocess, although 
as a check simulations including dd were run which confirmed expectations or lowering 
the asymmetry mildly. In this limit, the asymmetry is ApB — —0.5. 

First let us look at the events with ^/s = 24 TeV. We apply a cut rapidity 
Y > Yc = 1.0. In this way the misidentification due to the methodology is relatively 
small (~ 0.1). Also, the number of events that pass the cut is large enough as to make 
the statistical element very small (~ 0.03). As a result, the simulated set plotted in 
Fig 7 clearly shows a large asymmetry, compatible with the true one, and ~ Aa away 
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Mass (TeV) 

Figure 5. The discovery reach at a machine energy of =14 TeV and an integrated 
luminosity of 100 fb~^ for the specific model in the text. Orange (shaded) regions of 
the plot represent areas of the specific model that are either already excluded, where 
Drell-Yan production is suppressed or not accessible by the specific model. The (thick) 
blue line represents the minimum coupling ratio R required for at least 10 events to 
be observed. (Thin) blue lines follow the contours of an upper and lower Ict PDF 
uncertainty . The region between the upper and lower la uncertainties is shaded blue. 
The (thick) red line represents the minimum coupling ratio R for an asymmetry to be 
visible at the Str level (see text for definition). (Thin) red lines follow the contours of 
an upper and lower la PDF uncertainty (the lower line largely falls below the axis). 
The region between the upper and lower la uncertainties is shaded red. 

from zero. 

By comparison, in Fig. 8 we show a set of simulated events with the same choices 
of parameters, but a/s = 14 TeV. In this case we apply a lower cut of F > Fc = 0.75. As 
a result, the asymmetry of the simulated data is Aps = —0.45 ± 0.07(stat) ± 0.14(sys), 
lower than the expected true value, and affected by large uncertainties. 

A final comment is in order. In this numerical analysis we did not take into account 
the effect of experimental errors, nor any kinematical cuts that might be applied to 
extract analyzable data. For example, when | cos^| ~ 1 the experimental extraction of 
the kinematics of the two final leptons might be problematic. In order to understand 
how important this effect is, we extracted the asymmetry from the simulated data by 
imposing the further constraint | cos 6*1 < 0.8. This procedure effectively reduces the 
forward-backward asymmetry to ApB = —0.37 {ApB = —0.28) in the 24 (14) TeV case. 

These two examples show how the method works in practice. However caution 
must be noted that in performing such simulations there is a degree of sensitivity to 
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Mass (TeV) 

Figure 6. The discovery reach, at a machine energy of y/s =24 TeV and an integrated 
luminosity of 100 fb~^ for the specific model in the text. Orange (shaded) regions of 
the plot represent areas of the specific model that are either already excluded, where 
Drell-Yan production is suppressed or not accessible by the specific model. The (thick) 
blue line represents the minimum coupling ratio R required for at least 10 events to 
be observed. (Thin) blue lines follow the contours of an upper and lower Ict PDF 
uncertainty . The region between the upper and lower la uncertainties is shaded blue. 
The (thick) red line represents the minimum coupling ratio R for an asymmetry to be 
visible at the Str level (see text for definition). (Thin) red lines follow the contours of 
an upper and lower la PDF uncertainty (the lower line largely falls below the axis). 
The region between the upper and lower la uncertainties is shaded red. 

the choice of PDF sets. We illustrate here that the method functions but very different 
results may be obtained depending on how the large-^ distributions are modelled, and 
ultimately this problem will be resolved only once actual data are collected. 

6. Non-degenerate Couplings 

The region at low rapidity exhibits an interesting feature. Due to the differences in the 
up and down PDF distributions, the rapidity distributions of these two types of partonic 
events have very distinct shapes: up-type events peak at finite Y, while down-type peak 
at y = 0. In principle, this might allow to discriminate between different values of the 
couplings R-y' and Rz' of the new resonances to the up and down quarks. 

To investigate the feasibility of such a measurement we examine the two extreme 
cases, that Ry = versus Rz' = 0. We normalize the distributions to the same area, 
hence fixing the event number. The result is shown in Fig. 9. We see the largest effect is 
around Y = where there is a 1 part in 25 difference, or 5%. The chosen situation is the 
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Figure 7. An example simulation of the angular distribution for the predicted number 
of events at a yfs —2A TeV machine with L — 100 fb~^, for boson mass m = 2.6 
TeV and coupling ratio R — 0.4. The simulated data has been restricted to rapidity 
F > Yc = 1.0 and the angular distribution in the parton-parton centre of mass 
frame prediction based on the detector angular distribution is plotted. The observed 
asymmetry is Afb = —0.51 ± 0.03(stat) ± O.l(sys). This is the remnant of initial 
asymmetry of -0.42 that the data was simulated with. 



most extreme, a light mass in terms of what the precision constraints allow and a high 
energy collider. We can conclude that the effect of a non-degenerate R will be so small 
as not to be measurable, even assuming perfect detector efficiency and no kinematic 
cuts. 

7. Conclusions 

We performed a model-independent study of the Drell-Yan processes at the LHC for two 
new neutral gauge bosons with the quantum numbers of the photon and Z boson. We 
assumed the behavior of the new states are analogous to their standard-model cousins, 
but with couplings to SM fermions rescaled by a common factor of i? < 1. We also 
chose to look at the case where the peaks of the two new states were so close as to be 
indistinguishable in the invariant mass plots. 

We assumed perfect detector performance, omitting experimental errors, such as 
energy resolution, and without including the effect of kinematical cuts. We did not 
include any QCD corrections to the tree-level partonic calculation of pp — )■ jX^n' and 
pp — 7- e^e~ . In relation to the forward-backward asymmetry, we included in the analysis 
the main sources of theoretical uncertainty: a statistical element related to the number 
of expected events, a systematic due to the misidentification of the direction of the 
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Figure 8. An example simulation of the angular distribution for the predicted number 
of events at a -y/i =14 TeV machine with L — 100 fb~^, for boson mass m = 2.6 
TeV and coupling ratio R = 0.4. The simulated data has been restricted to rapidity 
Y > Yc = 0.75 and the angular distribution in the parton-parton centre of mass 
frame prediction based on the detector angular distribution is plotted. The observed 
asymmetry is ApB = —0.45 ± 0.07(stat) ± 0.14(sys). This is the remnant of initial 
asymmetry of -0.42 that the data was simulated with. 
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Figure 9. The expected rapidity distribution of a sample of events in the case that 
a single new Z' (blue) and a single 7' (red) is the only exchange particle present in a 
new resonance. The new boson has a mass set at 2.6 TeV and a machine energy of 
yfs = 24 TeV is used. The same number of events is assumed for both. 
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initial (anti-) quark, and the uncertainty due to the PDF distributions. 

The results can be interpreted as a best case scenario at the LHC, limited only by 
theoretical considerations. We compared different choices of the integrated luminosity 
and center-of-mass energy of the pp collisions. For illustration purposes, we also referred 
to one specific model, in which such mass-degeneracy of the new weakly-coupled gauge 
bosons is required by imposing the bounds on precision electro- weak parameters. 

For a 10 TeV collider, and luminosities in the 2 — 40 fb~^, bounds on new resonances 
no higher than 3 TeV can be obtained. With a 24 TeV collider the region well above 4 
TeV can be explored, depending on the specific coupling R. The uncertainty imparted 
upon this prediction from the PDFs is quite modest for a high energy collider, while 
at lower energies, like 10 TeV, the PDF uncertainty at large momentum fraction is an 
important limiting factor. This is in substantial agreement with the literature. 

We studied in detail the forward-backward asymmetry as a method of telling a 
two particle resonance from a one particle peak. The collider ability to measure an 
asymmetry at the 3cr level is heavily restricted to hght masses. Even at very high 
luminosities measuring the asymmetry for masses of 3 TeV is not possible at a 14 
TeV machine. Factoring in the PDF uncertainty, it is hard to conceive of any possible 
asymmetry measurements in a 14 TeV (or lower) collider, in relation to the scenarios 
considered here. By contrast, at a 24 TeV machine the PDF uncertainty is markedly 
lower and the reach far larger. At such a machine it may be possible to measure 
asymmetries and greatly restrict currently open parameter space. 

In many phcnomcnological models one would wish to measure precisely the 
physics of heavy gauge bosons. Measuring the forward-backward asymmetry in events 
with neutral intermediate states provides a theoretically clean signal allowing one to 
discriminate the presence of more than one boson in a heavy resonance. Such a 
measurement is very difficult in practice from an experimental viewpoint. On top of 
the practicalities there are the large PDF uncertainties at large momentum fraction. 
Together with the scaling of the event rate with s, the situation suggests upgrading the 
LHC to a higher energy (besides increasing the luminosity) to learn about such new 
bosons, if the new bosons are heavy as is suggested by the precision tests. 

Acknowledgments 

We would like to thank S. de Curtis and M. Fabbrichesi for useful discussions. The work 
of MR is supported by the STFC Doctoral Training Crant ST/F00706X/1, and MP is 
supported in part by the Wales Institute of Mathematical and Computational Sciences 
and also by the STFC Grant ST/G000506/1. 

References 

[1] See for instance T. Appolquist, B. A. Dobrescu and A. R. Hopper, Phys. Rev. D 68, 035012 (2003) 
[arXiv:hep-ph/0212073], and references therein. 



Mass-degenerate Heavy Vector Mesons at Hadron Colliders 



22 



[2] See for instance N. Arkani-Hamed, A. G. Cohen, E. Katz and A. E. Nelson, JHEP 0207, 034 (2002) 
[arXiv:hep-ph/0206021]; 1. Low, W. Skiba and D. Tucker-Smith, Phys. Rev. D 66, 072001 (2002) 
[arXiv:hep-ph/0207243]. 

[3] S. Weinberg, Phys. Rev. D 19, 1277 (1979); L. Susskind, Phys. Rev. D 20, 2619 (1979); 
S. Weinberg, Phys. Rev. D 13, 974 (1976). 

[4] W. A. Bardeen, C. N. Leung and S. T. Love, Phys. Rev. Lett. 56, 1230 (1986); M. Bando, 
K. i. Matumoto and K. Yamawaki, Phys. Lett. B 178, 308 (1986); B. Holdom and J. Terning, 
Phys. Lett. B 187, 357 (1987), Phys. Lett. B 200, 338 (1988); W. D. Goldberger, B. Grinstein 
and W. Skiba, Phys. Rev. Lett. 100, 111802 (2008) [arXiv:0708.1463 [hep-ph]]. 

[5] M. E. Peskin and T. Takeuchi, Phys. Rev. D 46, 381 (1992). 

[6] R. Barbieri, A. Pomarol, R. Rattazzi and A. Strumia, Nucl. Phys. B 703, 127 (2004) [arXivihep- 

ph/0405040]. 

[7] J. Hirn and V. Sanz, Phys. Rev. Lett. 97, 121803 (2006) [arXiv:hep-ph/0606086], JHEP 
0703, 100 (2007) [arXiv:hcp-ph/0612239]; D. K. Hong and H. U. Yee, Phys. Rev. D 74, 
015011 (2006) [arXiv:hep-ph/0602177]; M. Piai, arXiv:hep-ph/0608241, arXiv:hep-ph/0609104, 
arXiv:0704.2205 [hep-ph]; C. D. Caronc, J. Erhch and J. A. Tan, arXiv:hep-ph/0612242; 
M Fabbrichesi, M. Piai, L. Vecchi arXiv:0804.0124 [hep-ph]; K. Haba, S. Matsuzaki and 
K. Yamawaki, arXiv:0804.3668 [hep-ph]; J. Hirn, A. Martin and V. Sanz, arXiv:0807.2465 [hep- 
ph]: D. D. Dietrich and C. Kouvaris, arXiv:0809.1324 [hep-ph]. 
[8] For bounds in the context of the RandaU-Sundrum models, see for instance C. Csaki, J. Erlich and 
J. Terning, Phys. Rev. D 66, 064021 (2002) [arXiv:hep-ph/0203034]; K. Agashe, R. Contino, 
L. Da Rold and A. Pomarol, Phys. Lett. B 641, 62 (2006) [arXiv:hep-ph/0605341]; M. S. Carena, 
E. Ponton, J. Santiago and C. E. M. Wagner, Nucl. Phys. B 759, 202 (2006) [arXiv:hcp- 
ph/0607106]; Phys. Rev. D 76, 035006 (2007) [arXiv:hep-ph/0701055]; C. Bouchart and 
G. Moreau, Nucl. Phys. B 810, 66 (2009) [arXiv:0807.4461 [hep-ph]]; K. Agashe et al, Phys. 
Rev. D 76, 115015 (2007) [arXiv:0709.0007 [hep-ph]]. 
[9] "ATLAS detector and physics performance. Technical design report. Vol. 2," 

[10] M. Dittmar, Phys. Rev. D 55, 161 (1997) [arXiv:hep-ex/9606002]. 

[11] E. Accomando, S. De Gurtis, D. Dominici and L. FedeU, arXiv:0807.5051 [hep-ph]. 

[12] A. Martin and V. Sanz, arXiv:0907.3931 [hep-ph]. 

[13] M. Piai, arXiv:hep-ph/0608241. 

[14] M. R. Whalley, D. Bourilkov and R. G. Group, arXiv:hep-ph/0508110. 

[15] W. T. Giele, S. A. Keller and D. A. Kosower, arXiv:hep-ph/01G4052. 

[16] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. M. Nadolsky and W. K. Tung, JHEP 0207, 

012 (2002) [arXiv:hep-ph/0201195]. 
[17] M. Piai, arXiv:0704.2205 [hep-ph]; 



